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Abstract. 

Generation of terahertz electromagnetic radiation due to coalescence of upper- 
hybrid waves in the long- wavelength region of strong plasma turbulence driven by 
a high-current relativistic electron beam in a magnetized plasma is investigated. 
The width of frequency spectrum as well as angular characteristics of this radiation 
for various values of plasma density and turbulence energy are calculated using 
the simple theoretical model adequately describing beam-plasma experiments at 
mirror traps. It is shown that the power density of electromagnetic emission at the 
second harmonic of plasma frequency in the terahertz range for these laboratory 
experiments can reach the level of 1 MW/cm 3 with 1% conversion efficiency of 
beam energy losses to electromagnetic emission. 



PACS numbers: 52.25.Os, 52.35.Ra, 52.50.Gj 
1. Introduction 

It is well known that a turbulent plasma is a source of electromagnetic emission 
at the fundamental plasma frequency u> p and its second harmonic. These emission 
mechanisms have been found to play an important role in different space phenomena 
such as type III solar radio bursts [TJ [21 [31 IH [5] and radiations of planet's 
magnetospheres [BJ |7j- Plasma emissions at oj p and 2oj p have been also observed 
in laboratory beam-plasma experiments [H [H [101 E] • In contrast to the problem of 
type III radio bursts, in which the density of wave energy W is saturated at weakly 
turbulent levels W/nT ~ 10~ 5 (n is the plasma density, T is the temperature of 
plasma electrons), we focus our attention on the regime of strong plasma turbulence 
W/nT ~ 10~ 2 — 1CP 1 , which is more appropriate for laboratory experiments with 
powerful electron beams O EH E] . 

Beam-plasma experiments at mirror traps [141 115j have demonstrated that the 
relativistic electron beam with the typical energy 1 MeV and current density 1 kA/cm 
injected into the plasma with n = 2 ■ 10 14 cm -3 loses about 30% of its energy over the 
length of 1 m. It means that the averaged power density pumping by the beam to the 
plasma turbulence can be estimated as ~ 10 MW/cm 3 . The rate of beam energy losses, 
however, is not uniform over this length. The estimate for the peak pumping power 
can be derived assuming that a significant part of beam energy goes to excitation of the 
large amplitude coherent wave-packet at the entrance to the plasma column. Since the 
beam is decelerated over the packet length I ~ v^/T by Aw ~ v b T/uj p [TB], where T is 
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the growth-rate of the two-stream instability and Vb is the beam velocity, the maximum 
power of beam energy losses reaches the value ~ 100MW/cm 3 . Thus, if only 1% of 
pumping power is able to convert to electromagnetic radiation, the maximum power 
density of plasma emission produced in beam-plasma experiments at mirror traps can 
reach the range of 1 MW/cni . The same level of conversion has been observed recently 
in laboratory experiments [IT] with the low plasma density n = 2 ■ 10 14 cm~ 3 , that is 
why it allows us to suppose that, in a denser plasma, a beam excited turbulence can 
be also considered as an efficient source of powerful electromagnetic radiation. The 
aim of this paper is to calculate the conversion efficiency as well as absolute values 
of emission power for the case when second harmonic plasma emission falls in the 
terahertz frequency range. 

Our calculations of the spectral power of electromagnetic emission produced in 
a turbulent magnetized plasma are based on the theoretical model [17] . in which it 
is assumed that electromagnetic waves are generated predominantly in the source 
region of strong plasma turbulence due to coalescence of upper-hybrid waves. Since 
theoretical predictions have been found to agree with recent experimental results 
obtained at the GOL-3 multimirror trap in the low density regime, we can use this 
model to study spectral and angular characteristics of electromagnetic emission for 
the whole ranges of plasma density and turbulence energy, which can be achieved in 
our beam-plasma experiments. 

Thus, in Section 2, we formulate the main ideas of theoretical model that takes 
into account not only spontaneous generation of electromagnetic waves in coalescence 
processes, but also contributions of induced inverse processes. In Section 3, we present 
calculations of second harmonic emission power in the terahertz frequency range and 
study whether this emission can escape from the plasma. Our main results are 
summarized in concluding Section 4. 



2. Theoretical model 



To calculate the emission power from turbulent magnetized plasma we use the model 
of strong plasma turbulence proposed in [T7]. According to this model, we assume 
that most of wave energy is concentrated in upper-hybrid modes, which occupy the 
long- wavelength region of turbulent spectrum with wavenumbers k < yW '/nT ' /td 
(r_o is the Debye length). The saturation level of this energy is determined by the 
balance between the constant power P b pumping to the turbulence by the beam and 
the power dissipated due to the wave collapse. It results in the following relation 
between the pumping power and the turbulence energy: 

V2 / w \ 2 

where rrij is the ion mass and m e is the electron mass. 

We calculate the electromagnetic emission power P and the conversion efficiency 
e = P/Pb assuming that electromagnetic waves are generated due to coalescence of 
long- wavelength upper-hybrid waves. Since in the source region of turbulent spectrum 
these waves are not trapped in collapsing caverns, nonlinear interaction between 
electromagnetic and upper-hybrid modes can be treated in the framework of weak 
turbulence. In this case, generation of electromagnetic radiation is described by the 
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equation 



— ri. 



(2) 



dt ~ K ~ IK k ' 

where P k is the power of spontaneous emission produced in the coalescence process 
£ + £ — > t and 7fc is the nonlinear dissipation rate due to the decay t — > £ + £. 
In dimensionless units uj p t, uj/uj p , xlo p /c, kc/uj p , for time, frequency, position and 
wavenumber, respectively, the spectral wave energy is normalized by the condition 

W° k e?k = ^ 3 

nm e r 

where c is the speed of light. In the cold plasma limit, dimensionless values of P k and 
7fc take the forms 
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where we introduce the following notations: 
A*(k,u,) = |k.e£| 2 -fc 2 + u; 2 (e* k °i%e° k ) , 
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Here, eigenfrequencies wjl and eigenvectors e£ of linear plasma modes are determined 
by the standard dispersion equation with the dielectric tensor i k (I is the unit matrix) . 
In contrast to the similar calculations of second harmonic emission [T51 IIH] based on 
the standard weak turbulence theory, we take into account model damping of two- 
time correlation functions, which is used to describe the effect of finite life-time of 
upper-hybrid plasmons due to their scattering off density fluctuations with the typical 
frequency v = u) p W e / (nT). It results in correlation broadening of the resonance 
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0, which is described in (|4j) by the function 
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This function shows that the width of frequency spectrum of second harmonic 
electromagnetic emission depends essentially on nonlinear effects as well as on thermal 
and magnetic corrections to the linear dispersion of upper-hybrid modes. In order 
to minimize the frequency width of radiation, we should consider the case, when 
the ratio of the electron cyclotron frequency to the plasma frequency becomes low 
Q = uj c /ujp — 0.2 resulting in comparable contributions of magnetic and thermal 
effects for the typical plasma temperature T = 1 keV in mirror traps. As to account 
for the effect of finite temperature, we modify the eigenfrequencies and eigenvectors of 
linear plasma modes, but we neglect modifications in the nonlinear current G\ u k k , 
which is much less sensitive to thermal corrections than Ak,ki,k 2 - Thus, to calculate 
w£, we use the fluid approximation and solve the linear dispersion relation with the 
dielectric tensor: 

/ k 2 V 2 \ 

(11) 

(12) 

(13) 
(14) 

Sy Z — —E Z y — —I A - , (15) 

J k 2 V 2 + tt 2 \ 
e zz = l-Ail- J, (16) 

A = (uj 2 - fl 2 - k 2 VZ + %kf,VS 

\ LO 2 11 

where Vj, = 3T/(m e c 2 ), magnetic field is directed along z-axis and k = (k±, 0, fcy) is 
the wave vector with the length k = [kj_ + k 2 ) 1 / 2 . 

If the path lengths of spontaneously generated electromagnetic waves are larger 
than the typical size of confined plasma Ik — v g /jk > L (v g is the group velocity of 
electromagnetic wave) , decay processes t — > £ + 1 do not play a role, and the second 
term in ([2| can be omitted. Thus, in the case of azimuthally symmetric turbulence, 
the spectral emission power in units of nm e c 2 is given by the integral 

7T 

dP f ( k 2 \ 

1 ' 

where k(oj) is the solution of u> = u)\ and 9 is the polar angle of k. 




3. Computation results 

Let us study the spectral and angular characteristics of second harmonic plasma 
emission for various regimes of beam-plasma interaction, which can be realized in 
mirror traps. We will vary the plasma density from 2 • 10 14 cm' 3 to 5 • 10 15 cm~ 3 for 
different fixed values of turbulence energy W/nT = 0.01, 0.05, 0.1 and for the fixed 
parameters T = 1 keV and ft = 0.2. 
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Figure 1. The spectral power dP/ (dwd6) (in W/cm 3 ) for x-mode emission (a) 
and o-mode emission (b). The path length (in cm) for x-modes (c) and o-modcs 
(d). 

We assume here that the source region of turbulent spectrum contains the 
anisotropic population of resonant waves, which are directly pumped by the beam, 
and the isotropic population of nonresonant waves, which are the products of beam- 
driven waves scattering off density fluctuations. In our computations, the isotropic 
part occupies uniformly the spectral region kc/uj p £ (0.1, k m c/ui p ), where the upper 
bound corresponds to the typical wavenumber of the modulation instability k m ~ 
y/W/nT ' /rjj and the lower bound excludes modes with wavelengths larger than the 
typical plasma size. According to the model [T7], the anisotropic population of beam- 
driven waves contains a small fraction of energy (10%) and occupies the region: 
kc/uj p e (1.1, 1.3) and 6 € (0, 0.3). 

In the regime with n — 3 • 10 15 cm~ 3 and W/nT — 0.05, computation results for 
the spectral density of emission power dP / (dujdO) for extraordinary (x) and ordinary 
(o) electromagnetic modes are presented in Figures [l|a) and |T|(b) , respectively. It is 
seen that x-mode emission is dominated by the contribution of obliquely propagating 
waves with 9 = 30° and 9 = 150°, whereas the most intensive o-mode emission 
is concentrated near the transverse to the magnetic field direction. Calculations 
of the path length Ik for these emissions are shown in Figures [ijc) andgd). As 
one can see, the path length reaches the minimal value 20-40 cm for longitudinally 
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propagating modes regardless of their polarization. For transverse propagation, this 
length increases up to 90 cm for the o-mode and exceeds 3 m for the x-mode. It means 
that emissions of both electromagnetic modes, generated in the plasma column with 
the diameter 5-6 cm typical to our beam-plasma experiments, are able to escape from 
the plasma and can be used for different applications. 

Let us now find out how the frequency spectrum of electromagnetic emission, 
described by dP/duj, depends on the plasma density. As a function of dimensionless 
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Figure 2. The spectral power of second harmonic emission dP/dui in MW/cm 3 
for different values of turbulence energy W/nT. 
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Figure 3. Dependences of / (a) and A/ (b) on the plasma density for different 
values of turbulence energy W/nT. 

frequency the form of this spectrum is affected by the level of turbulence only. 

For the density n — 3 • 10 15 cm" 3 and different values of turbulence energy, examples 
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of frequency spectra, derived by summation over polarizations, are shown in Figure 
[2] In dimensional form, the linear frequency / corresponding to the maximum of 
the spectral power dP/du and the width at half maximum A/ for a fixed turbulence 
energy depend on the plasma density as follows: /, Af oc n 1 / 2 . Figures [3|a) and|3|b) 
demonstrate that / and especially A/ grow with the increase of the turbulence level. 
It means that the spectral width of plasma emission in regimes of rather strong plasma 
turbulence is mainly determined by correlation broadening of resonances in nonlinear 
three-wave interactions. 

To compute the integral power of second harmonic emission from the unit volume 
of turbulent plasma, we integrate the function dP/dui over the frequency range 
uj/uj p £ (1.6, 2.6) and sum over polarizations. It is seen from Figure|4j that this power 



1000 , 



6 100 



10 r 



0,1 



co /co =0.2 
c p 

T =1 keV 



-W/nT=0.1 
- W/nT=0.05 
W/nT=0.01 



10 

,„H -3 

n, 10 cm 



Figure 4. The dependence of integral second harmonic emission power on the 
plasma density for different values of turbulence energy W/nT. 

increases like P oc n 3 / 2 with the increase of plasma density for the fixed turbulence 
energy and reaches the value 1 MW/cm 3 in the regime with n = 3 • 10 15 cm -3 and 
W/nT = 0.05. It is interesting to note that these calculations adequately describes 
recent experimental results obtained in the strong magnetic field fl — 0.8 [TT]. Indeed, 
the 1 kW/cm 3 range of emission power, observed experimentally in the regime with 
n = 2 • 10 14 cm -3 and W/nT — 0.01, is also reproduced by calculations presented in 
Figure [2] From comparison with the results of [T7] we can conclude that the external 
magnetic field broadens the spectral width of electromagnetic emission, but does not 
change drastically the integral power. Since the pumping power ([!]) has the same 
dependence on plasma density as the emission power, the conversion efficiency e in 
our model is completely determined by the turbulence energy. For W/nT = 0.05, this 
efficiency is estimated as e ~ 1%. 

To make clearer understanding how computation results are sensitive to the 
shape of turbulent spectrum, we calculate the spectral emission power for different 
spectra corresponding to the same turbulence energy (Figure [5]) . When we replace 
the uniform isotropic spectrum of upper-hybrid modes used in previous calculations 
on more realistic falling spectra, changes in the integral emission power do not exceed 
20%. 
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Figure 5. Different shapes of the isotropic part of turbulent fe-spectra (a) and 
the corresponding spectral power of electromagnetic emission (b). 



It should be emphasized that the high value of turbulence energy W/nT = 0.05, 
which is required to generate 1 THz emission with the power density 1 MW/cm 3 , 
can be achieved in beam-plasma experiments at the GOL-3 multimirror trap. Indeed, 
for the typical beam energy 1 MeV and current density 1.5 kA/cm 2 , the pumping 
power reaches the required level of 100 MW/cm inside the region of most intensive 
beam-plasma interaction, where large amplitude coherent wave-packets are excited. 
Since the length of this region is estimated as I ~ Vb/T ~ 1 — 3 cm, the total power of 
terahertz emission in our laboratory experiments can reach 30-100 MW. 



4. Conclusion 



The spectral power of second harmonic electromagnetic emission, generated in the 
long- wavelength region of strong plasma turbulence due to coalescence of upper-hybrid 
waves, was calculated for different regimes of beam-plasma interaction, which can be 
realized in modern experiments at mirror traps. Computation results have shown 
that beam energy transferred to plasma oscillations in the regime of strong plasma 
turbulence can be converted efficiently to electromagnetic radiation. We have studied 
spectral characteristics of plasma emission in wide ranges of plasma density and 
turbulence energy and shown that the power of terahertz radiation in our beam-plasma 
experiments can reach tens of MW. Taking into account the angular distribution of 
this emission and possibility to change the emission frequency by varying the plasma 
density we come to the conclusion that terahertz radiation generated in a beam-driven 
strongly turbulent plasma can be very attractive for different applications. 



Generation of powerful terahertz emission 



9 



Acknowledgment s 

Authors thank LA. Kotelnikov for useful discussions. This work is supported by 
grant 11.G34.31.0033 of the Russian Federation Government, President grant NSh- 
5118.2012.2, Russian Ministry of Education and Science and RFBR grants 11-02- 
00563, 11-01-00249. 

References 

[1] Gurnett D A, Anderson R R 1976 Science 194 1159 

[2] Kruchina E N, Sagdeev R Z, Shapiro V D 1980 JETP Lett. 32 419 

[3] Robinson P A, Cairns I H and Willes A J 1994 Astrophys. J. 422 870 

[4] Mel'nik V N and Kontar E P 2003 Sol. Phys. 215 335 

[5] Li B, Willes A J, Robinson P A and Cairns I H 2005 Phys. Plasmas 12 012103 

[6] Gurnett D A, Shawhan S D and Shaw R R 1983 J. Geophys. Res. 88 329 

[7] Cairns I H and Menietti J D 2001 J. Geophys. Res. 106, 29515 

[8] Benford G, Tzach D, Kato K, Smith D F 1980 Phys. Rev. Lett. 45 1182 

[9] Hopman H J and Janssen G C A M 1984 Phys. Rev. Lett. 52 1613 

[10] Baranga A B, Benford G, Tzach D, Kato K 1985 Phys. Rev. Lett. 54 1377 

[11] Arzhannikov A V, Burdakov A V, Kuznetsov S A, Makarov M A, Mcklcr K I, Postupaev V V, 
Rovenskikh A F, Sinitsky S L, Sklyarov V F 2011 Fusion Sci. and Technol. 59 (IT) 74 

[12] Vyacheslavov L N, Burmasov V S, Kandaurov I V, Kruglyakov E P, Mcshkov O I and Sanin A 
L 1995 Phys. Plasmas 2 2224 

[13] Vyacheslavov L N, Burmasov V S, Kandaurov I V, Kruglyakov E P, Meshkov O I, Popov S S 
and Sanin A L 2002 Plasma Phys. Control. Fusion 44 B279 

[14] Arzhannikov A V et al. 2003 JETP Lett. 77 358 

[15] Arzhannikov A V, Burdakov A V, Koidan V S, Vyacheslavov L N 1982 Physica Scripta T2, 303 

[16] Timofeev I V, Lotov K V 2006 Phys. Plasmas 13 062312 

[17] Timofeev I V 2012 Phys. Plasmas 19, 044501 

[18] Willes A J and Melrose D B 1997 Sol. Phys. 171 393 

[19] Kuznetsov A A 2007 Plasma Phys. Rep. 33 482 



